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Abstract 

Delicate cancellations among diagrams can result in the vanishing of threshold am- 
plitudes in the standard model. This phenomenon is investigated for multi-Higgs 
production by scalar, vector, and fermionic fields. A surprising gap is found where 
although the tree-level amplitude for 2 incoming gauge bosons to produce N v particle 
may vanish, the amplitude for the production of N v + 1 particles is nonzero. 



One of the most interesting pieces of phenomenology to come from the study of 
multi-particle processes is amplitude nullification!!! @ . The interest in this phe- 
nomenon was spurred by an observation about the tree-level amplitude for two on- 
shell scalar particles scattering to produce n particles. At the kinematical threshold, 
this 2 — * n amplitude vanishes for all n > 4 in A0 4 theory with an unbroken reflection 
symmetry©. In A</> 4 theory with broken symmetry, it was noted that the tree-level 
threshold amplitude of 2 — > n production vanishes for all n > 20. 

The nullification phenomenon forbids the threshold amplitude for one off-shell 
scalar decaying into n scalars, 1 — ► n, from developing an imaginary part at the 
one-loop level. At the one-loop level, any relevant unitary cut must cross exactly two 
lines. Any such cut will relate the imaginary part of the amplitude to the product of 
the vanishing 2 — > n amplitude. 

At the threshold, amplitudes of 1 — ► n processes still exhibit factorial growth, 
and, therefore, appear to violate unitarity. It is not clear what happens above the 
threshold, and it may be possible that diagram cancellation plays a critical role in the 
restoration of unitarity. It can be speculated that this nullification is the manifestation 
of some as to yet undiscovered symmetry. There exists a scenario where nullification 
could provide a clue as to the restoration unitarity at high multiplicity in the standard 
model©. In order for this to occur, the coupling of the Higgs field to the other fields in 
the standard model would have to take on very special values. Requiring these special 
couplings to occur would place restrictions on the ratios of masses in the standard 
model, and be an attractive way of explaining the particle mass spectrum. 

For simplicity, we will consider a U{1) gauge theory coupled to a fermionic field 
and a complex scalar field. All of the features in the following discussion can be 
easily generalized to the standard model. The masses of the gauge and fermionic 
fields are generated in the standard way through the Higgs mechanism. The mass 
and quadratic coupling of the physical scalar particle shall be denoted as m and A/4, 
respectively. After symmetry breaking, the scalar field develops a vacuum expectation 
value v = \Jtri 2 / (2A). For simplicity, we will work in units where m = 1. The 
amplitude calculations shall rely on the reduction formula techniqueB for calculating 
multi-particle amplitudes. The key observation for using this technique, is that at 
the threshold, all fields have zero momenta and are functions of time only. The LSZ 
reduction simplifies to a one dimensional problem. Using this method, the amplitude 
for the threshold decay of one scalar to n scalars is, 



1 



( 



d_ 

dz 



) 



n 



(n\a\0) 



(0+ |a(x)|0-) 



(1) 



where z(t) is the response of the free Higgs field in the presence of a source, p, and is 
proportional to e imt , and (0 + | cr| — ) is the response of the interacting Higgs field in 
the presence of a source. The tree-level amplitudes are generated by approximating 
the response by the solution to the classical field equations, cr c i(t). This well known 
generating function for 1 — > n processes is0, 



The first loop correction to the amplitude is generated by expanding the mean 
value of the field around the classical solution, (0 + |<j(x)|0— ) = cr c i(x) + (Xi(x), where 
(Ji(x) is the mean value of the quantum part of the field. The equations of motion 
can be expanded around <r c i(x), and only the leading contribution in A needs to be 
retained at the one-loop level, 



where the three functions, V'(t), represent the contributions to the amplitude due 
to scalar, fermion, and vector field loops. Each of the contributing terms, V(t), in 
equation ([3D contains a term representing the equal time two-point field correlation 
function evaluated in the a c i background field. This can be visualized as the evaluation 
of a tadpole Feynman diagram, with the equal time Green function coming from the 
propagator in the loop (for a discussion on using diagrammatic techniques in the 
calculation of multiparticle processes, see ||). 

This two point Green function can be found by taking the inverse of the second 
variation of the Lagrangian evaluated with o = a c i. The Green function is a transfor- 
mation of the amplitude of 2 incoming particles to produce n scalars at the threshold. 
For physical processes, the interesting object is the double pole of the amplitude as 
the momenta of the incoming particles are taken to be on mass shell. These processes 
can be characterized by a single parameter, u = \/k 2 + 1 = n/2, where k is the spa- 
tial momentum of the incoming particles, and n is the number of particles produced 
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in the final state. If the Green function does not have a pole for u = n/2, then the 
tree-level amplitude for 2 — > n particles is zero. 

In the following sections, we will discuss how particles from each sector of the 
standard model contribute to the one-loop amplitude of 1 Higgs — ► n Higgs and 
where nullification in the production of Higgs particles occurs in this sector. 

The Higgs Sector 

Both the one-loop correction to 1 — > n processes and 2 — > n nullfication for 
incoming scalars have been discussed previously in [£|. The first variation of the 
scalar potential contributes terms to the equation of motion that are both quadratic 
and cubic in the scalar field, a. Expanding around the mean field, 




VLalarit) = 3\/ £ Mx)<T q (x)) + 3\* d (<T q (x)<T q (x)), (4) 

where {a q (x)a q (x)) is the spatially independent two-point Green function in the clas- 
sical background field, G(x,x') = (Ta q (x)a q (x')), taken in the limit x — > x' . As 
discussed in a previous work®, the Green function is the inverse of the second vari- 
ation of the Lagrangian evaluated in the classical background field. 




G(x,x') =-i5(T,r'). (5) 



The component of G(x,x') with spatial momentum k, denoted as G u (t,t j ) has 
been constructed from the homogeneous equation with the operator in equation (||) 
B. For future reference, we express G u (t, t') in terms of a general Green function, 



G>,t,t') = l/(2uj)(f^T)f-(T')Q(r-r') + f-(r)f^T)Q(r'-r)), (6) 
where r is defined by, 

e T = ~fy(t), (7) 

and the two solutions to the homogenous differential equation with the operator in 
equation (|5|) are 

ft = e^ T F(-a, a + 1; 1 ± 2u; -^ T ), (8) 

1 + e T 
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and F(a, b; c; z) is the hypergeometric function. The hypergeometric series in equation 
(0) will terminate when a is an integer. For incoming scalar particles, G w (t,t j ) = 
G w (2,t, r'). If the series did not terminate, the Green function would have poles 
for all half integer values of u. Since the series terminates, there are poles only for 
u) = ±1/2, ±1. Thus, the tree-level amplitude for two on-shell scalars to form n 
scalars vanishes at the threshold for n > 2. 

Summing over all momenta modes, the equal time Green function is, aside from 
a renormalization of m and A, 



<T q \X t 



aq(x)) = . (9) 

Mi - *(*) A^) 4 

The Higgs sector cannot contribute an imaginary part to the one-loop 1 — > n 
amplitude. There is no two body phase space when uj = 1/2, 1 so it is not possible 
to make a nonzero unitary cut to contribute to the imaginary part. 

Fermionic Contribution 

The fermionic contribution to the equation of motion is found by varying the 
fermionic portion of the Lagrangian with respect to a. As in the case of the scalar 
field, the mean value of a is expanded around the classical solution, a c i(t). The mean 
value of all other fields is set to zero. The fermionic term in equation @ is, 

v' fermion {t) = ^(*M*)>, (io) 

where (ijj(x)ijj(x)) is the trace two-point fermionic Green function, (Ti()(x)i(j(x')) , and 
m f is the mass of the fermion field as generated by the Higgs mechanism. The Green 
function has been discussed in and can be found by calculating the inverse of the 
variation of the Lagrangian with respect to ip{x) and ip(x) evaluated in the classical 
background field, 



Vflt 

m f a d (t) 



(TiP(x)il;(x')) = -i5\x - x'). (11) 



The spinor structure of equation (|Tl|) is simplified by the introduction of the Green 
function H(x, x') defined by, 
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and analytic continuation into imaginary time as in the previous section. Since spatial 
momentum is conserved, the Green function, H( be constructed from the 

solution to the homogeneous differential equation, 

d 2 n / o m 
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/(r) = (13) 

The matrix 7 in equation ([IB] ) has eigenvalues of +1 and —1. As was noted by 
Voloshin in 0, in order for nullification to take place, the coefficient of the sech 2 (r/2) 
term must be N(N + l)/4 with integer N for both eigenvalues. This imposes the 
conditions on m/, 



4m 2 f + 2m f = N f (N f + l), (14) 



and, 



4mJ - 2m f = N' f (N' f + 1), (15) 

for some integer values of Nf and iVj. Both equation ([14]) and equation (|15D are 
satisfied when 2m f = Nf and Nf = Nf — 1. 

In a basis where the matrix H is diagonal, two of the diagonal elements of H will 
be G w (Nf, t, t') and two will be G w (Nf — 1, r, r'). The fermionic Green function for 
a general fermionic coupling is, 
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Equation ([TB]) contains two types of terms. One type of terms contributes only to 
the renormalization of m and A. The other terms actually contribute to the physical 
correction to the 1 — > n amplitude. In the theory with Nf = 1 (mj = m^/2), only 
the contribution to the renormalization of the coupling constant and mass is present. 

The theory with the simplest nontrivial (^}(x)ip(x)) and diagram nullification is 



characterized by iV 



2, i.e. 



m* 



mft. In this theory, the contributions from 



Guj(Nf — l,r, r') contribute only to the renormalization of m^ and A. Apart from 
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renormalization terms, the fermionic correlator arising from G u {2, t,t') is, 

ty(x)${x)) = ? JU ■ 17 

3tt 2 (1 -z(t)J\/2) 5 

It is worth noting that the fermionic contribution to the equations of motion and, 
hence, to the one- loop correction is exactly — | that of the scalar contribution. 
Gauge Field Sector 

The situation with vector fields is more complicated. The qualitative behavior of 
2 vector — > n scalar processes is very dependant on the polarization of the incoming 
gauge particles. Care must be taken in the choice of gauge. Nonphysical amplitudes 
like those of the 1 — > n processes discussed in previous sections are gauge depen- 
dent. While it is easiest to calculate in unitary gauge, the gauge loop corrections 
to processes like 1 — > n are divergent. These divergences will cancel when a gauge 
independent physical amplitude is calculated. We will work in unitary gauge for 
simplicity of calculation, but we will concentrate on discussing the properties of the 
physical amplitude of two incoming vector bosons scattering to form many scalars. 

The two point Green function is calculated as in previous sections by taking the 
inverse of the second variation of the Lagrangian. In the unitary gauge, the two point 
vector Green function, (V p (x)V v (x')) , satisfies the condition, 



d 2 + ^(a d (t)+vr 

v 



(V^x)V„(x')} - d»d<>(V p (x)V u (x')) = -ig^5\x - x'). 



There are no additional ghost contributions in unitary gauge. It is convenient to 
rotate to a frame of reference where the spatial momentum is along the " 1" axis. The 
spatially fourier transformed Green function, 

GV(M') = / (V,(x)V u (x'))^d\ (19) 
has only two nonvanishing derivitives, do and d\. For the transverse components, G22 



and G33, the last term on the left hand side of equation (18) does not contribute. The 
Green functions may be written as in previous sections in terms 

of the solutions to the homogeneous differential equation constructed from the 
operator in equation (p~8f) , 



GW = W (j£{T)f- (t')Q(t - t') + /- (r)/+ (r')e(r - r')), (20) 
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Where W is the Wronskian of the two solutions. In analogy with previous sections, 



we define the variable, N v , which satisfies, 



4m: 



For the the transverse components, 



N V (N V + 1). 



(21) 



m = m = e" T F(-N v , N v + l,lT 2u; — -). (22) 

The Wronskian for these solutions is 2u. The transverse part of the Green function 
will have poles at all half integer values of to unless N v is an integer. As noted by 
VoloshinB, the physical tree-level amplitude for 2 incoming transversely polarized 
vector particles to produce n scalars at the threshold vanishes for all n > N v . The 
transverse off-diagonal elements of G^ u are zero. 

For longitudinally polarized particles, things are a little bit more complicated. 
Differentiation of equation (|18|) yields the auxiliary equation, 



0, 



M 2 tanh 2 r/2 



G^ = d v 5(r-r'). 



Searching for solutions of the form, 

1 
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(23) 



(24) 



tanhrtanhr'" v ' ' ' M£ tanh 2 r / '2 
the Green function, H(t,t'), becomes the object of interest. Combining equations 



and (|i~8D with the definition of H(t, t'), 



d 



uo 2 -m 2 tanh 2 (T/2) 



2sinh r/2 



H(r,r') 



{to 2 - iy 



5(t-t'). (25) 



As in equation (|20|), H(r) can be constructed from the two solutions to the ho- 
mogenous differential equation with the operator in equation (|25|). The operator is 
the same as the one that appears in the Poschl- Teller equation. The solution as 
described in the Appendix is, 



± 



-2— + atanh(r/2) + coth(r/2) 

UT 



e^F{-{N v -l),N v ;l±2u; 



1 + e T 



). (26) 
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Note that the normalization of the Green function, H(t,t'), on the right hand 
side of equation ( pop is not unity. As with the case of transverse gauge particles, there 



are poles in /# at every half integer value of uj when N v is not an integer. When N v 
is an integer, there are no poles for u > (N v — l)/2. The Wronskian of the functions 
fij can be conveniently evaluated at the point r = oo. The result is 



Au 2 - (N v + 1) 2 J . (27) 

There is an additional pole contributed by the Wronskian at the point to = (N v + 
l)/2. The unusual pattern of the nullification is the following: all process are allowed 
for n < N v . The production of n = N v + 1 (which is forbidden in the transverse 
mode) is allowed. There is a small gap in the allowed region, when N v particles are 
produced, which is not allowed. Prouction of n particles for n > N v + 1 is forbidden. 
The remaining components of G^ v can be easily constructed from the above formulae, 
and have the same poles as Goo- 

Because of the pole at n — N v + 1, gauge field loops can contribute an imaginary 
part to 1 — > n amplitude. The amplitude for the 1 — > n process is not a gauge 
invariant object, and is highly divergent in the unitary gauge. 

Concluding Remarks 

The phenomenology associated with nullification can be quite complex. It's pres- 
ence can prevent the one- loop 1 — > n amplitude from developing an imaginary part. 
It can also cause the cancellation of 2 — * n diagrams in all sectors of the standard 
model for n > N for some value of N. 

The behavior in the longitudinal vectors is quite unexpected. While for incoming 
transverse vectors there is cancellation of 2 — > n diagrams for all n > N v , the process 
with n = N v + 1 is allowed for longitudinal polarized particles. In addition, the 
process with n = N v is forbidden for longitudinal polarized particles. The pole at 
n = N v + 1 is what allows the 1 — > n amplitude to develop an imaginary part. It is 
of interest to note that the number of poles is the same regardless of the polarization 
of the incoming particles. 

Although there exists much knowledge about the phenomenology of nullification, 
little is known about its root cause. It is possible that it is associated with some 
symmetry in the standard model that may play a part in the restoration of unitarity 
at large multiplicities. However, other scenerios have been suggested!!. 

I would like to think Mikhail Voloshin for many inciteful discussions. This work 
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Appendix: Solution of the second Poschl- Teller Equation through the 
methods of supersymmetric quantum mechanics 

In the section on vector fields, the Green function was related to the solution to 
the Poschl- Teller equation, which has the form, 



d 2 a(a + 1) 
dx 2 cosh 2 fx) 



/ = A/. 



(2f 



sinh (x) 

For the Green functions of vector fields, 6 has the value of 1. 

This problem may be elegantly solved by using supersymmetric arguments (many 
reviews have been written on the uses of supersymmetry in quantum mechanics, e.g. 
H). Consider supersymmetric quantum mechanics with two supercharges defined to 
be, 

Q± = [±^ + W{x)\<f k (29) 



where a denotes the appropriate 2x2 Pauli spin matrix and the choice of potential is 
made to be, 



W(x) = atanh(x) + 6coth(x) (30) 
The supersymmetric Hamiltonian is the anticommutator of the two supercharges, 



H = -{Q + ,Q_} 



d 2 



dx' 



+ (a + b) 2 ^— (a - a 3 ) + , 2 (6 - a 3 ). (31) 

cosh (x) sinn (x) 



This operator <j\ had two eigenvalues, +1 and —1. These degenerate eigenstates 
shall be denote as if>±, and obey the eigenvalue equations, 



d 2 



+ (a + b) 2 - 



a(a =f 1) 6(6 + 1) 



(32) 



dx 2 v cosh z (x) sinh z (x) 

Consider the case 6=1. The second term in H + disappears to leave a familiar 
operator. The solutions are known to be, 

1 



ij,+(x) = e ±UJX F(-(a -1),o;1tw; 



1 + e a 



(33) 



where u — J E — (a + l) 2 , and F(a, b; c; z) is the hypergeometric function. By virtue 
of the superymmetry, the eigenvalues of H~ are simply Q_%p + . It is easy to see the 
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the solution to equation (28) for b = 1 is simply 
/ 



— : — h a tanh (x) + coth (x) 

ax 



±v^ j p ( _ (a _ 1);a . lTv ^.^ :) _ ( g 4) 

1 + e x 



Solutions to equation (|2"8"D for larger integer values of b may be contructed by iterating 
this procedure. 
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